On the spatial structure of anomalously localised states in disordered conductors 
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We study the spatial structure of wave functions with exceptionally high local amplitudes in the 
Anderson model of localisation. By means of exact diagonalisations of finite systems, we obtain 
and analyse images of these wave functions: we compare the spatial structure of such anomalously 
localised states in quasi-one-dimensional samples to that in three-dimensional samples. In both cases 
the average wave-function intensity exhibits a very narrow peak. The background intensity, however, 
is found to be very different in these two cases: in three dimensions, it is constant, independent of 
the distance to the localisation centre (as expected for extended states). In quasi-one dimensional 
samples, on the other hand, it is redistributed towards the localisation centre and approaches a 
characteristic form predicted in [A. D. Mirlin, Phys. Rep. 326, 249 (2000)]. 



Statistical properties of physical observables in disor- 
dered electronic quantum systems have attracted consid- 
erable interest in the last decade. In such systems, quan- 
tum interference may cause the (non-interacting) conduc- 
tion electrons to localised. In three dimensions (3D), this 
Anderson localisation occurs when the disorder strength 
exceeds a critical value. Beyond this value (which de- 
pends on the Fermi energy and the symmetries of the 
system), electron wave functions are confined to a lim- 
ited spatial region of the sample. 

In the metallic regime, on the contrary, wave func- 
tions typically spread over the whole sample, and they 
contribute to electron transport. However, some wave 
functions show localised behaviour even in this weakly 
disordered regime (these are so-called anomalously lo- 
calised states, abbreviated as ALS in the following). This 
leads to a non-zero, albeit small, probability of observ- 
ing exceptionally large wave-function amplitudes, often 
in the form of a log-normal distribution function of wave- 
function intensities. ALS in electronic conductors have 
been studied intensively in recent years, using the. 
called diffusive non-linear sigma model (DNLSM)o 
An overview of the main results and predictions based 
on the DNLSM is given in Ref. §. Moreover, possible 
complications due to non-diffusive, so-called ballistic ef- 
fects on length scales smaller than the mean free path 
are discussed. As was pointed out in Ref. ||, these may 
modify the predictions of the DNLSM (see also Refs. ^ 
and |Io| ) . ALS are expected to occur in lower-dimensional 
disordered systems, too, when the disorder is weak. 

These interesting analytical results have motivated 
a number of numerical studies: in Ref. 0, for exam- 
ple, log-normal statistics of wave-function amplitudes in 
two-dimensional (2D) conductors near the delocalisation- 
localisation transition was observed. In Ref. [n] it was 
confirmed that as the disorder is reduced to reach the 
weakly disordered regime, the distribution function re- 
mains log-normal. It has, however, not beenjjpssible to 
resolve a discrepancy (between the DNLSMtm and the 
so-called direct optimal fluctuation methodta) in the pre- 



diction of the parameters of this distribution. As far 
as the numerical work is concerned, the situation may 
be summarised as follows. While it can be concluded 
that the DNLSM appears to be adequate in the quasi- 
one-dimensimal (Q1D) Anderson model under certain 
conditionald'El, the DNLSM rna y nat correctly describe 
ALS in two and three dimensionsEjO, at least for the pa- 
rameter values considered in_these studies. On the other 
hand, it was found recentlyEj that the DNLSM appears 
to describe the statistics of rare events adequately in a 
2D kicked rotor. 

A reason for the possible failure of the DNLSM to de- 
scribe ALS in the 3D Andersosumodel may be the im- 
portance of short length scales™, see also Ref. [l0| The 
DNLSM is based on the semiclassical picture of a diffus- 
ing electron in a random potential. In this picture the 
smallest relevant length scale is the electronic mean free 
path I between elastic collisions. Therefore the DNLSM 
cannot describe situations where length scales smaller 
than t are important. Such a situation could occur if 
ALS were created not through semiclassical diffusion,djujt 
through local potential wells trapping the electw»sBli3. 
The mechanism giving rise to ALS is expectedoEni3 to 
crucially depend on the dimensionality of the system and 
may determine their spatial structure. 

Finally, within the DNLSM, it is possible to obtain the 
wave- function statistics directly in Q1D, using a transfer- 
matrix technique. In 2D and 3D, on the other hand, a 
further saddle-point approximation is necessary. 

Given this situation, further numerical work describing 
ALS in disordered conductors is called for. In the follow- 
ing, we describe results of exact diagonalisations of finite 
Anderson models of localisation, yielding averages of the 
spatial structure of ALS. While this spatiaLstructure has 
been studied in detail within the DNLSMErlij, numerical 
imaging of wave functions with anomalously amplitudes 
has not yet been performedEl 

The Anderson model is defined by the tight-binding 
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Hamiltonian 



H 



,t„,cl c, + 



(1) 



on a hypercubic lattice. Here cj and c r are the creation 
and annihilation operators at site r, the hopping am- 
plitudes are \t TT '\ = 1 for nearest-neighbour sites and 
zero otherwise. The on-site potentials v r are Gaussian 
distributed with zero mean and (iytv) — (W 2 /12)5 rr ' . 
As usual, the parameter W characterises the strength of 
the disorder and (• • • ) denotes the disorder average. We 
study finite Q1D and 3D lattices. 

It has been suggested in Ref . || to characterise the spa- 
tial structure of ALS by means of conditional averages of 
the form (V q \tp(r)\ 2q ') t . Here (•••){ denotes an average 
over all wave functions with t — V|^>(0)| 2 , q = 1,2,..., 
and V is the volume of the system (wave functions are 
normalised so that (\ipjM\ 2 ) = V~ 1 ). For large values of 
t, the average (y q \ip(r)\ 2q ) describes how wave functions 
decay, on average, away from the localisation centre. 

In the metallic regime, typical wave fun ct i on s fluctu- 
ate as described by random matrix theorjOlla (RMT). 
Depending on the symmetries of ([!]), Dyson's Gaussian 
orthogonal or unitary ensembles are appropriateta. We 
refer to these two cases by assigning, as usual, the pa- 
rameter j3 = 1 to the former and /3 = 2 to the latter. 
Within RMT, for r f 0, 

' 1 if q = 1 

(V q \ip{r)\ 2q ) t = { 3 if? = 2andj9 = l l (2) 
2 if q = 2 and = 2 . 

This reflects that in a metallic system, wave functions 
spread uniformly over th£_ whole sample with spatially 
short-ranged correlationsctu. In the presence of ALS, the 
conditional averages (y q \ip(r)\ 2q ') are expected to differ 
from (§). 

In order to characterise the spatial structure of ALS in 
the weakly disordered Anderson model, the conditional 
averages (• ■ ■ )t were calculated within the DNLSM in 
Refs. H and |l6[ For q = 1 the authors write 

(V\^{r)\ 2 ) t =g{E,t;r)/f{E,t), (3) 

where f(E,t) is defined as 

f(E lt ) = A(Y,5(t-V\^(0)\ 2 )) (4) 

3 

Here A is the mean energy level spacing and (• • • ) denotes 
a combined disorder and energy average (over a small 
interval of width n centered around E). The function 
g(E, t; r) is defined as 



g(E,t;r) = A(j2V\Tp j (v)\ 2 6(t-V\i> j (0) 



Ei~E 



(5) 



Below, we take the ^-functions in (0) and (|^) to be 
slightly broadened with small but finite 7 > 0. 



In close vicinity of the localisation centre, ALS are 
found to exhibit a very narrow peak (of width less than £). 
The expressions below apply for r > £ and thus describe 
the sppoth background intensity, but not the sharp g@aJc 
itselM For (3 = 2 one obtains for a Q1D conductor!!^ 



f(E,t) = ^ [W^(t/X,T + )W^(t/X,T^ 
and (assuming |r| = r > £) 



g(E,t;r) = -X- 



WW(t/X, Tl ,T 2 )WW(t/X,T_) 



f 



(6) 



(7) 



where t + = [L — x)/£, r_ = x/£, t± = r/£ and T2 = 
(L — x — r)/£. Here, X = L/£, L is the length of the 
sample, and £ is the localisation length. Moreover, x is 
the distance of the observation point from the edge of the 
sample (c.f. Refs. [|f[6|). The function W^(z,t) obeys 
the differential equation 



OT 



(8) 



with initial condition W^'(z, 0) = 1. The function 
W^(z,t, t') obeys the same differential equation, but 
with the initial condition W (2) (^,0, r) = zW (1) (z,t). 

For large values of t, it is suggested in Ref. |l^ that the 
background intensity should be given by 



(9) 



where in accordance with the above r > £ is assumed, and 
also r <C £■ Corresponding expressions may be obtained 
for higher values of q. One expectsEJ (for £ < r <C £) 



(10) 



In 3D, b«-|Contrast, results corresponding to (8-10) are 
not knowntj, but it is expected that the background in- 
tensity of ALS is characteristically different froaubhat in 
Q1D samples: ALS in 3D systems are expectedtm to ex- 
hibit a very narrow maximum near the localisation centre 
(the width of this peak is a matter of Cttcsent debateEJH). 
The background, however, is expectedtlli3 to decay very 
quickly to (V\ip(r)\ 2 ) t ~ 1. Furthermore, fluctuations 
around this average are expected!^ to be described by 
RMT [see Eq. @}. 

In the following we-,describe and discuss results of 
exact-diagonalisationsEJ of the Anderson Hamiltonian (|l|) 
and compare to the results of the previous section. We 
emphasise that at least in 3D it is important to use a 
Gaussian distribution for the on-site potentials in order 
to be able to compare to the analytical predictions: ALS 
are possibly non-universal, their spatial structure may 
depend on the properties of the random potential. 
Our numerical results for jY\ip(j)\ 2 ) t and (V^ 2 |^(r) | 



are summarised in Figures [lj, |^, and [| . Figs, [l] and 3 
show results for Q1D samples with (3 — 2, while figs. 2 
and H show results for 3D samples and = 1- 
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FIG. 1: Spatial structure of anomalously localised wave func- 
tions with t — \ip(0)\ 2 V. Solid lines: Numerical results for 
V = 128 x 4 x 4 (Q1D case), disorder W = 1.6 and en- 
ergy E ~ —1.7, averaged over 40 000 wave functions. Dashed 
lines: Analytical predictions with X — 0.97 [full formula, with 
Eqs. (^) and (fit)]. The dash-dotted line shows the asymptotic 
formula, Eq. (0), for t = 25. 




FIG. 2: Spatial structure of anomalously localised wave func- 
tions with t = \if)(0)\ 2 V. Numerical results for V = 48x48x48 
(3D case), disorder strength W — 2 and energy E ~ — 1.7, av- 
eraged over 88 wave functions. The symbols show the numer- 
ically calculated values for t = 1 (o), t = 10 (□), £ = 20 (+), 
and t = 30 (*). The line shows the constant RMT average 
intensity (\ip(r)\ 2 V) t = 1. 



In all cases we see a very narrow peak at the locali- 
sation centre. Our main results, however, concern char- 
acteristic differences between the distribution of back- 
ground intensities in Q1D and 3D samples. In the Q1D 
case we observe a global redistribution of background in- 
tensity towards the localisation centre (as compared with 
uniformly spread RMT wave functions). The numerical 
results are very well described by Eqs. (6-8). The asymp- 



FIG. 3: Solid lines: The function (V 2 \ip(r)\ A ) t in the Q1D 
case averaged over the same wave functions as in Fig. 
The dashed lines show the function 2[(|V"-i/;(r)| 2 ) t ] 2 , where 
(\Vip(r)\ 2% j t was calculated from eqs. (8-10). 



totic formula (9) considerably underestimates (V\ip(r)\ 2 ) t 
for the values of t used in Fig. fij. In the 3D case, by 
contrast, the background intensity is roughly constant 
(Fig. ||). This is consistent with the qualitative picture 
summarised above. 

Figs. |||andU show the second moments, i.e., the case 
q — 2. In Q1D samples, Eq. ( |Io| ) appears to be valid for 
large t, as expected. In 3D samples, the fluctuations of 
the background intensity appear to be consistent with the 
RMT statistics (although the scatter is large). That is, 
they are roughly constant as a function of r, as predicted 
m Ref. [l| 

However, a closer inspection shows (inset of 
Fig. ||) that the second moment is somewhat higher than 
expected according to RMT. Thkis possibly due to the 
finite conductance in the systenoEj. The dip at r ~ 7 
with t = 30 could be due to insufficient averaging. Our 
current data do not permit to draw a definite conclusion 
here, it could also represent a systematic effect. 

Summarising our results for r > I, we have made the 
following observations. First, in our simulations, the 
spatial structure of ALS depends on the dimension. In 
Q1D samples, we have observed a global redistribution of 
the background intensity towards the localisation centre 
(as compared with typical, uniformly spread extended 
states). In 3D samples, on the other hand, the back- 
ground intensity is roughly constant; the wave-function 
intensity is significantly increased only in the very narrow 
vicinity of the localisation centre. 

Second, the DNLSM appears to describe the spatial 
structure of background intensity of ALS adequately in 
the Q1D Anderson model. The agreement is good for all 
values of t studied here. This verifies the assumption that 
the semiclassical picture of a diffusive electron accounts 
for the origin of ALS in the Q1D Anderson model. 



4 



100 m 




FIG. 4: The function (V 2 |i/>(r)| 4 ) t in the 3D case. The sym- 
bols show the numerically calculated values. The parameters 
and symbols are the same as in Fig. El The line shows the 
constant RMT fluctuations, (V 2 |i/>(r)|*) t = 3. 



Third, we must emphasise that at least in the Q1D case 
it appears to be difficult to reach the asymptotic regime 
where (||) is valid. We hasae thus not been able to observe 
the line shape suggested^ to be characteristic of ALS in 
Q1D samples. We expect, however, that this character- 
istic shape is approached further as t is increased. 

Fourth, fluctuations around the average ALS, in Q1D, 
are consistent with Eq. (|l0|). The agreement is the bet- 
ter the larger £ is, as seen in Fig . |^. This is so because 
Eq. ( |l0| ) was derived (in Ref. jig ) using asymptotic ex- 
pressions. 

Fifth, in 3D samples the fluctuations around the av- 



erage ALS appear to be consistent with RMT, although 
the second moment is somewhat larger than expected ac- 
cording to RMT. 

Finally, we briefly consider the region r < £. In the 
3D case, the central peak is found to be very narrow. 
It is possibly much narrower than the electronic mean 
free path I which is of the pcder of several lattice spac- 
ings in the metallic regimeOj. This may indicate that 
the DNLSM is not applicable in this case, since £ is the 
shortest relevant length scale in the DNLSM. It may also 
point towards the conjecture of Refs. || and ^| where high 
wave-function amplitudes were suggested to arise from 
partial trapping of electrons in rare local potential cavi- 
ties. Such features are not included in the DNLSM. At 
this point, however, we cannot offer quantitative results 
concerning the parametric dependence of the width of 
the central peak. 

In conclusion, we have studied the spatial structure 
and statistics of anomalously localised, states in the An- 
derson model. Our results indicateO that the spatial 
structure of ALS in Q1D and 3D samples is very dif- 
ferent, as surmised in Refs. p||l6|, and ^. Our results are 
consistent with the idea that the origin of ALS is different 
in Q1D and in 3D. In order to decide to which extent lo- 
cal potential traps are relevant in 3D, it would be of great 
interest to study the parametric dependence of the width 
of the local maximum seen in Fig. |^. In this context, a 
continuous model would probably be more suitable than 
the discrete lattice Hamiltonian ([!]) studied in the present 
paper. 
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